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Context. The absolute calibration of space-borne instruments in the soft X-ray regime rests strongly on model spectra of hot white dwarfs. 
Aims. We analyze the Chandra LETG+HRC observations of the white dwarfs HZ43 A and Sirius B and of the neutron star RX J185635-3754 
with the aim of resolving current uncertainties in the soft X-ray spectral fluxes and photospheric parameters of the three stars. We apply the 
derived photon spectra to a cross-calibration of the LETG+HRC-S with the short-wavelength EUVE spectrometer and the ROSAT PSPC. 
--^ Methods. We tie HZ43 A to the flux of RXJ1856 in the 44-48A range and perform a simultaneous least squares fit to the LETG+HRC 
i-C ' spectra of the three stars. This allows us to determine an internally consistent set of spectral energy distributions and an empirically derived 
wavelength-dependent correction to the LETG+HRC-S effective area. We employ NLTE model atmospheres calculated with TMAP for the 
Q ■ white dwarfs and a two-blackbody model for RX J 1856, tied to the respective optical fluxes. 

S-H Results. The two-blackbody model for RX J1856 features a hot spot on a cooler star and yields kr spol = 62.8 ± 0.4 eV and kT stM = 32.3 ± 
(/} , 0.7 eV with a stellar angular radius as seen from infinity of 0.1371 ± 0.0010km pc~'. For HZ43A, our fit yields 7^=51 126+660 K and 

log g=7. 90+0.080 (cgs) with anti-correlated errors (1-cr) which include not only the statistical but also the systematic uncertainties of the fit. 
^ HZ43AB displays a previously detected bremsstrahlung component with a temperature kT ^ 0.6 keV. For Sirius B, we find T^g =24923+1 15 K 
• i-H , for fixed logg=8.6. The calibration of the short- wavelength EUVE spectrometer differs from that of the LETG+HRC-S by 15 + 7%. The ROSAT 

PSPC is found to be correctly calibrated within a few percent and reports of a major miscalibration are unfounded. 
jH Conclusions. We have obtained improved parameters for RX J185635-3754, HZ43 A, and SiriusB which fit the observations from the optical 
. . i to the soft X-ray regime. Our approach allows us to quote their absolute spectral fluxes at selected wavelengths which may aid the calibration 

of other space-borne instruments. 

Key words. Methods: data analysis - Stars: white dwarfs - Stars: neutron - Stars: fundamental parameters - Stars: individual (HZ43A, 
SiriusB, RX J185635-3754) - X-rays: individual (HZ43 A, SiriusB, RX J185635-3754) 



1. Introduction 

Hot white dwarfs with pure hydrogen atmospheres are well 
established standard stars in the optical and ultraviolet (e.g. 
Bohlin et al. 2001) and have played a decisive role in the cali- 
bration of satellite instruments in all spectral regimes from the 
optical to the extreme ultraviolet and soft X-rays. The FOS 
and STIS on the Hubble Space Telescope (Bohlin 2000), the 
Extreme Ultraviolet Explorer (EUVE) (Boyd et al. 1994), and 
the Low Energy Transmission Grating Spectrometer on board 
Chandra (LETG+HRC-S) with the High Resolution Camera 
as detector (Brinkman et al. 2000; Pease et al. 2000) are promi- 
nent examples. The latter two instruments were calibrated on 



* Based on observations collected with the Chandra Observatory, 
the Hubble Space Telescope, the EUVE satellite, and ROSAT. 



the ground, but their calibration was partly relinquished in 
favor of an in-flight calibration using hot white dwarfs. The 
Position Sensitive Proportional Counter (PSPC) on board the 
Rontgen Satellite (ROSAT), on the other hand, was accurately 
calibrated on the ground (Triimper 1982; Pfeffermann et al. 
1987) and not recalibrated in flight. Based on observed PSPC 
spectra of white dwarfs, however, Napiwotzki et al. (1993); 
Jordan et al. (1994); Wolff et al. (1995, 1996) claimed the 
PCPC calibration at soft X-ray energies to be incorrect by a 
factor of about two. Given this important role of white dwarfs 
it is deplorable that their soft X-ray fluxes are much less well 
documented than their optical and ultraviolet fluxes. Given, fur- 
thermore, the large body of scientific data of all types of objects 
potentially affected, the resolution of such open calibration is- 
sues is of general interest. 
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Our initial aim was to investigate the PSPC calibration 
which requires, however, to take recourse to spectrally higher 
resolved data as provided by the EUVE or the Chandra 
LETG+HRC spectrometers. Since the effective area of the lat- 
ter is itself uncertain by as much as 25% at intermediate wave- 
lengths around 60A (e.g. Pease et al. 2003), we found it nec- 
essary to perform a more detailed cross-calibration of the three 
instruments. As a result, we are able to determine accurate soft 
X-ray fluxes of the hot white dwarfs HZ43 A and Sirius B and 
of the neutron star RX J 185635 -3754 (henceforth RX J 1856) 
and thereby establish them as standards in the soft X-ray 
regime 1 . Although one of the two white dwarfs (HZ43A) is 
a well known standard for wavelengths A > 912A (Bohlin et 
al. 2001) and the optical flux of the other (SiriusB) is reason- 
ably well measured (Barstow et al. 2005), the extrapolation of 
the optically adjusted model spectra into the soft X-ray regime 
suffers from the wagging-tail problem. We solve this problem 
by tying the soft X-ray tail of the spectrum of HZ43 A to that 
of RX J 1 856. This adds a second fix to the spectrum of HZ43 A 
besides that in the optical and effectively constrains its spectral 
shape. 

The organization of the paper is as follows. In Sect. 2, we 
describe the analysis of the LETGS and PSPC data and list 
the EUVE archive data used. Sect. 3 is devoted to the descrip- 
tion of our theoretical white dwarf spectra and summarizes the 
available information on T e ff and log g of HZ43 A. Our prin- 
cipal results are provided in Sect. 4. It starts with information 
on the optical fluxes to which our model spectra are normal- 
ized (Sect. 4.1). In Sect. 4.2, we derive an accurate spectrum of 
RXJ1856 that serves as the short-wavelength fix for HZ43 A. 
In this context, we discuss interstellar absorption near the car- 
bon K-edge at 43. 7 A to some detail. The essence of the paper 
is then presented in Sect. 4.4, the simultaneous fit of RX J 1 856, 
HZ43A, and SiriusB that results in an internally consistent 
set of photon spectra and simultaneously in the wavelength- 
dependent correction functions to the LETG+HRC-S effective 
areas at wavelengths A > 44A. Here, we discuss also the sys- 
tematic errors of our approach. The stellar parameters implied 
by our fit are presented and discussed in Sect. 4.5. Sect. 4.6 
compares our LETGS calibration to that of the EUVE short 
wavelength spectrometer and Sect. 4.7 finally demonstrates the 
accurate calibration of the ROSAT PSPC. The conclusions are 
summarized in Sect. 5. 



2. Observations and Data Analysis 

We have extracted all observations of HZ43A, SiriusB, and 
RXJ1856 taken with the Chandra LETG+HRC-S, the EUVE 
spectrometer, and the ROSAT PSPC from the respective 
archives. Table 1 contains a log of the observations. 



2.1. Chandra LETG+HRC spectra 

The spectra were extracted from the reprocessed level 2.0 
event files. The corresponding total exposure times are listed 
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1 We consider neither quasars and BL Lacs because of their inherent 
variability nor supernova remnants because of their extended nature. 



Fig. 1. Long-term variation of the integrated level 2.0 LETG+HRC 
count rates of HZ43 A for two spectral intervals. The two sets of tick 
marks indicate the times when the LETG+HRC spectra of SiriusB 
and RX J 1856 were taken. 



in Table 1 . The extraction region recommended in the Chandra 
Proposers' Observation Guide was used and the background 
determined from wide regions below and above the source ex- 
traction area. Further details are given by Burwitz et al. (2003). 
Since the background is variable, we investigated the influence 
of rigorously discarding intervals of high background on the 
resulting count rate spectra. Our tentative elimination scheme 
(in short referred to as level 3) reduces the exposure times of 
RX J1856, HZ43 A, and Sirius B to 66%, 84%, and 61% of the 
level 2.0 accepted times, respectively, and is more restrictive 
than the condition which Drake et al. (2002) imposed on the ac- 
cepted telemetry rate. Nevertheless, the wavelength-integrated 
count rates increase by only 0.8%, 0.1%, and 6.8%, respec- 
tively. These values measure the additional dead-time resulting 
from the time intervals of enhanced background. We opt to ac- 
cept the level 2.0 spectra with the small dead-time corrections 
applied to them as a multiplicative factor. Within their statis- 
tical errors, the level 3 and the dead-time corrected level 2.0 
spectra are indistinguishable frome one another. 

The long series of spectra taken of HZ43 A between 1999 
and 2006 lends itself to an analysis of the long-term behavior 
of the HRC sensitivity. The zeroth order count rate of HZ43 A 
decreases slowly with time at a rate of 0.98 + 0.15% per year 
or with a time scale of about 100 years. This decrease seems 
to depend on wavelength, with a low of the decay rate of 0.5% 
per year near 110A, rates of 0.9 ± 0.2% for 70-80A, and a 
more uncertain 2.4±0.8% for 44^-9 A (Fig. 1). For Sirius B and 
RX J 1856, the time bases are much more restricted as indicated 
by the tick marks in Fig. 1 . The decay rates of their zeroth order 
count rates are ~ 0.7% and close to zero, respectively. The spa- 
tial and temporal gain changes of the HRC-S are a known effect 
and a study for times before JD2452900 has been presented by 
Pease & Drake (2003). We do not attempt to correct for these 
changes, but rather chose a pragmatic approach. We place a si- 
multaneous spectral fit on a secure basis by selecting a subset 
of the HZ43 A exposures with an epoch nearly equal to that of 



K. Beuermann et al.: Soft X-ray standard stars HZ43 A, Sirius B, and RX J185635-3754 



3 



Table 1. Journal of observations extracted from the Chandra, EUVE, 
and ROSAT archives for the present analysis. Only the total exposure 
times are given, except for the ROSAT PSPC where the experimental 
setups differ and the individual exposure times are given. 



Object 



Obs. Identif. 



Year 



Exp. (ks) 



a) Chandra LETG+HRC-Spectrometer 

HZ43A 00059,01011,01012, 1999,2001 254.8 

02584, 02585, 03676, 2002 

03677, 05042, 05044, 2003, 2004 

05957, 05959, 06473 2005, 2006 

SiriusB 01421,01452,01459 1999,2000 64.1 

RXJ1856 00113,03380,03381, 2000,2001 501.9 
03382, 03399 

b) EUVE-Spectrometer 

HZ43A 1904N, 1620N 1994,1995 130.0 

1049N, 0030N 1996, 1997 

SiriusB 0147N, 0039N 1993,1996 246.0 

c) ROSAT PSPC 

HZ43A rpl00308 1990 21.5 

rf200418 1991 21.6 

rpl41916/17 1992 6.9 

SiriusB rf200422 1991 3.2 

RXJ1856 rp200497 1992 6.3 



the RXJ1856 observations, i.e., the exposure-weighted mean 
Julian day JD 2452130. We use the HZ43 A exposures 00059 
through 03677 which comprise 157239 s of exposure with an 
epoch different from that of RX J 1856 by only 38 days. A sim- 
ilar selection is not possible for Sirius B which is centered 508 
days earlier. We do not correct for the possible 1 % effect. 

We start using the first-order effective areas of the 
LETG+HRC-S in the positive and negative dispersion direc- 
tions as recommended in the November 2004 release of the 
Chandra X-ray Center 2 and refer to them as the nominal areas 
A°. Throughout this paper, we use the efficiencies £k of orders 
k = 2 to 6 relative to the first order as given in the same release. 
The areas A are reliable to better than 15% shortwards of the 
carbon K-edge at A — 43. 7 A, but may need corrections possi- 
bly as large as 25% in the 60-80A regime (Pease et al. 2003). 
We distinguish, therefore, A° from A = aA°, where a(A) is 
a wavelength dependent adjustment factor which equals unity 
shortwards of the carbon K-edge and is allowed to deviate from 
unity longwards of it. Our correction function a derived below 
is strictly valid only for the epoch JD 2452130 and may differ 
for other epochs during the lifetime of the HRC at the percent 
level. 

The calibration of the LETG+HRC-S at long wavelengths 
is based, in part, on Sirius B and HZ43 A (Brinkman et al. 2000; 
Pease et al. 2000), with a certain preference given to SiriusB 
because of its better defined log g (Barstow et al. 2005; Holberg 
et al. 1998). The drawback is the lack of a useful signal from 
Sirius B at A < 60A, while HZ43 A can be detected at least 



2 http://cxc.harvard.edu/cal/Letg 
(updated November 2004). 
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down to A = 44A. Since we wish to connect the spectrum of 
our white dwarf calibrator to the harder spectrum of RX J 1856, 
we give preference to HZ43 A. This has the added advantage 
that RXJ1856 and HZ43A are both observed on-axis while 
the Sirius pointings are slightly off axis. The latter may affect 
the simultaneous fit if the HRC efficiency differs spatially (in 
addition to temporally) for a given wavelength. We neglect this 
possible effect. 

We bin the LETG+HRC spectra to 0.5A in wavelength 
which is an adequate compromise between spectral resolution 
and statistical significance of individual data points in the less 
well exposed wavelength regions. The detector gaps for the ob- 
servations of RX J 1856 and HZ43 A are conservatively defined 
as 49-58A and 58-69A for the negative and positive disper- 
sion directions, respectively. For SiriusB, the gaps overlap and 
statistically significant data of exist only for A > 61 A. 

2.2. EUVE spectra 

We obtained six night time EUVE spectra of HZ43 A and two 
of Sirius B. The archived photon spectra were used and we re- 
frained from a re-extraction of the data. The archived spectra 
are based on corrected effective areas as described by Boyd et 
al. (1994). Two spectra of HZ43 A (9302180059, 940325231 1) 
were discarded because of excessive noise and low flux, respec- 
tively. The remaining four are listed in Table 1 . They differ in 
their general flux levels by +5%, but no systematic time vari- 
ation is seen. We did not use the very faint and ill-determined 
EUVE spectrum of RX J1856 discussed by Pons et al. (2002). 

2.3. ROSAT PSPC spectra 

All ROSAT spectra of HZ43 A, SiriusB, and RXJ1856 taken 
with the Position Sensitive Proportional Counter (PSPC) were 
extracted from the ROSAT archive and re-analyzed. Since the 
three sources are very bright, we used a generous extraction 
radius of 5 arcmin radius which was only reduced for the de- 
termination of the high energy bremsstrahlung tail of HZ43 A 
(O'Dwyer et al. 2003). The background was taken, if possible, 
from a surrounding ring free of additional faint sources. We list 
in Table 1 the subset of observations discussed in detail here, 
one each from the 1990 Program Verification (PV) phase, the 
1991 Announcement of Opportunity 1 (AO-1), and the 1992 
AO-2. The PV observations were taken with PSPC-C and the 
AO-1 observations with PSPC-B before the reduction of the 
high voltage on October 11, 1991. The detector response ma- 
trix for both data sets is DRMPSPC-AOl (Briel et al. 1995). 
For the AO-2 observations after the voltage change , the appro- 
priate matrix is DRMPSPC. In order to allow a straightforward 
comparison with the previous analyses of Wolff et al. (1996, 
and references therein), we do not apply any adjustments to 
the matrices and refrain from corrections for the known (small) 
temporal and spatial gain shifts of the PSPC (Prieto et al. 1996; 
Snowden et al. 2001). For SiriusB and RXJ1856 only one 
PSPC observation each was performed, of which the former 
is in AO-1 and the latter in AO-2. The observations with identi- 
fication numbers starting with rp were taken in the open config- 
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Fig. 2. NLTE pure hydrogen spectra for Tgg =50700 K and logg=8.0, 
calculated with TMAP using two versions of the bound-free and free- 
free absorption coefficients, Seaton's approximation (dashed curve) 
and the full Karszas & Latter (1961) description (solid curve). Both 
models yield practically identical results at optical/ultraviolet wave- 
lengths, but differ by a factor of 1.33 at 44A (vertical dotted line). 
The insert shows an expanded view. Two further dotted lines mark 
A = 1300A and A = 5450A (see text). 

uration, the ones starting with rf include the boron filter which 
reduces the count rate for the very soft white dwarfs by about 
a factor of eight as a result of the added absorption and the 
narrower width of the response in photon energy. 

All integrated PSPC count rates quoted in Sect. 4.7 below 
refer to energy channels 1 1-60 (of 256). We note already here 
that our count rates are lower than those presented by Wolff 
et al. (1996) who seem to have included channels below the 
recommended lower limit of channel 1 1 . 

3. Model atmospheres and the case of HZ43 

In this Section, we describe the non-LTE model atmospheres 
used in our analysis and summarize previous results on the pho- 
tospheric parameters of HZ43. 

3.1. NLTE models for pure hydrogen atmospheres 

We employ TMAP, the Tubingen NLTE Model Atmosphere 
Package (Werner et al. 2003; Rauch & Deetjen 2003), in its 
most recent version for the calculation of plane-parallel, static, 
pure-hydrogen models. Hydrogen is represented by a model 
atom with the lowest levels up to principal quantum number 
n = 14 treated in NLTE. The higher levels are dissoluted 
rapidly towards the inner photosphere at the high g considered 
here, e.g., to less than 9% occupation probability at a column 
density of only 10~ 5 g cirT 2 for n — 14 and the r e jf = 27900 K 
and logg = 8.6 of Sirius B (Hummer & Mihalas 1988). Hence, 
this number of NLTE levels is sufficient. Two additional levels 
are treated in LTE. The NLTE levels are fully coupled by ra- 
diative line transitions. An approximate formula for the Stark 
broadening (Unsold 1968) is employed for the calculation of 
their absorption coefficients. Collisional transitions are consid- 
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T e ff (10 3 K) 

Fig. 3. Summary of r^-log g determinations for HZ43 A from fits to 
the Balmer lines (circles), the Lyman lines (triangles), and the EUVE 
continuum (lying crosses): 1 = Napiwotzki et al. (1993); 2 = Finley 
et al. (1997); 3 = Barstow et al. (2003); 4 = this work (STIS Balmer 
lines); 5 = Dupuis et al. (1998) (ORFEUS); 6 = Barstow et al. (2001) 
(HUT); 1 = Barstow et al. (2001) (ORFEUS); 8 = Barstow et al. (2003) 
(FUSE); 9 = this work (FUSE); A = Barstow et al. (1995) (EUVE); and 
B = Vennes & Dupuis (2002) (EUVE). The error ellipses indicate the 
results of the present fits to the Chandra soft X-ray continuum and 
refer to the 68%, 90% and 95% confidence levels. 



ered between all levels. The atmospheres are calculated within 
optical depths of log r = -8 . . . + 4. In order to model the flux 
at the Lyman threshold precisely within 900 - 1000 A, a narrow 
frequency grid was used with a spacing equivalent to 0. 1 A. 
For the synthetic spectra, we use the line-broadening tables of 
Lemke(1997). 

The short-wavelength continuum benefits from the use of 
//-cross sections from Sutherland (1998) 3 which replace the 
older fit formula from Mihalas (1967). More importantly, the 
continuum flux depends on the treatment of the bf- transitions. 
Using the full formalism of Karzas & Latter (1961) raises the 
flux at 44A by 33% compared with spectra calculated with the 
older Seaton formula and Gaunt factor (see Fig. 2). We find 
very good agreement of our results with the short-wavelength 
continuum of Barstow et al. (2003). Given the sensitivity of 
the soft X-ray fluxes to details of the calculation, we caution 
against the indiscriminate comparison of different versions of 
model spectra used over the decades in the calibration of vari- 
ous satellite instruments sensitive in the soft X-ray regime. 

Another aspect suggested to be important at soft X-ray en- 
ergies is the suppression of the flux by Compton rather than 
Thomson scattering. Although Madej (1998), claimed a size- 
able effect, Suleimanov et al. (2006) show that it is only at the 
1% level near the carbon K-edge in HZ43 A and Sirius B. Since 
the Compton effect deforms the model spectrum, we investigate 
the effect quantitatively below. 



3 http://astro.uni-tuebingen.de/~rauch/TMAP/UserGuide/ 
UserGuide.html 
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3.2. The case of HZ43A 

Figure 3 summarizes previous measurements of T e g and log g 
of HZ43 A supplemented by our present results. The individ- 
ual parameter combinations scatter over 6000 K and 0.4 dex, 
respectively, with no preference for any individual method to 
prefer a particular region in parameter space. Systematic dif- 
ferences in the deduced parameters from Lyman and Balmer 
lines were extensively studied by Barstow and collaborators, 
who show that they are still small near 50000 K (Good et al. 
2004). The data points in Fig. 3 cluster around r eff = 50500 K 
and logg = 7.95, but it is not certain that such a poll bears 
physical significance. The cause for the partly disparate data 
points is likely systematic in nature and may have its origin in 
observational problems, e.g., the presence of stray light in the 
line centers; in differences of the fitting procedures, e.g., the 
way a quasi-continuum is adjusted; in differences of models, 
e.g., between LTE and non-LTE or between the individual ver- 
sions of the line broadening theory; and in physical differences 
as the unrecognized presence of trace metals and their levita- 
tion (Barstow et al. 2003; Bohlin 2000; Chayer et al. 1995; 
Napiwotzki 1997; Schuh et al. 2002). Clearly, the extremely 
good counting statistics in the spectra of HZ43 A in all spectral 
ranges are not the only error source. The error ellipses in Fig. 3 
denote our result for HZ43 A based on the interpretation of the 
Chandra soft X-ray spectrum. This result takes systematic er- 
rors partly into account and is described below. 

4. Results 

Model fits to the individual Chandra LETG+HRC spectra us- 
ing the nominal (November 2004) first order effective areas 
and the nominal efficiencies of the higher orders relative to the 
first order yield marginally acceptable reduced xl = 1 -45 and 
2.04 for RX J 1856 and Sirius B, respectively, but an enormous 
xl = 32.40 for HZ43 A. Following suggestions by others (e.g. 
Braje & Romani 2002; Drake et al. 2002; Pease et al. 2000, 
2003), we interprete these deficiencies as caused by the remain- 
ing errors in A and search for a way to remove these defects. 
To this end, we include both hot white dwarfs in a simultane- 
ous fit together with the neutron star RX J 1856, which allows 
us to tie HZ43 A to the flux level of the latter and to the better 
defined effective areas at A < 40A. In what follows, we discuss 
the normalizations of the model spectra and the treatment of 
the interstellar absorption which is of particular importance for 
RX J 1856 and subsequently describe the combined fit. 

4.1. Visual and EUV fluxes 

We require that our incident spectra comply with the find- 
ings in other wavelength regions. For HZ43A, we adopt 
the visual magnitude V = 12.909 (Bohlin 2000), or equiv- 
alently adjust the model to the HST/STIS spectral flux at 
5450A, /5450 = 2.558 x lO-^ergcirrV'A- 1 . In order to 
appropriately fix the interstellar absorbing hydrogen column 
density, we adopt a mean flux in the 475^-95A region of 
/485 = 1.84 x 10 _12 ergcm _2 s _1 A _1 , which is the spectral 
flux of the long-wavelength EUVE spectrum used by Wolff et 



al. (1999) and slightly exceeds that of Dupuis et al. (1995). 
For Sirius B, we adopt the 4600A HST/STIS flux / 4600 = 
2.860 x 10~ 12 ergcm _2 s _1 A _1 taken from the carefully cali- 
brated G430L HST/STIS spectrum of Barstow et al. (2005) 
which puts us on his optical flux scale. The mean 475^-95A 
EUVE flux of Sirius B, ~ 6 x 10~ 15 erg cirrV A" 1 , has a huge 
error, but a flux of this order is consistent with the interstel- 
lar neutral hydrogen column density reported by Hebrard et al. 
(1999), N m = (6.5 + 2.0) x 10 17 cirT 2 . ForRXJ1856, we use 
the visual flux determined by van Kerkwijk & Kulkarni (2001), 
/ 5000 = 2.96 x 10~ 19 erg cm _2 s _1 A~' at 5000A, and account for 
extinction using A v = _/V H i/(1.87 x 10 21 x H cirT 2 ) (Clayton et 
al. 2003), where Nm is the fitted interstellar column density of 
atomic hydrogen in cirT 2 and x H is the ionization fraction of 
hydrogen. 

4.2. Interstellar photoelectric absorption 

The strength of the interstellar carbon K-shell absorption edge 
in the spectrum of RXJ1856 is central to our discussion and 
requires that we discuss this aspect to some detail. We note that 
the application of the nominal LETG+HRC effective areas A 
yields an unphysical carbon jump of the RX J 1856 spectrum in 
emission, which corresponds to a negative carbon abundance 
a(C). The appropriate treatment of carbon absorption allows us 
to fix the correction to the LETG+HRC effective area on the 
long-wavelength side of the carbon jump. This is the first step 
in our approach to an internally consistent set of photon spectra 
for our three stars. 

In treating interstellar photoelectric absorption, we largely 
follow Wilms et al. (2000). We use the cross-sections of hydro- 
genic ions as given by the analytic formula from Spitzer's text- 
book, the cross sections of Yan et al. (1998) for neutral helium, 
and the cross sections of Balucinska-Church & McCammon 
(1992) and Verner & Yakovlev (1995) as implemented in the 
most recent version of XSPEC 4 for the heavier elements. For 
the short lines of sight to the objects considered here, we ne- 
glect molecules. We use the protosolar abundances of Lodders 
(2003) as representative of cosmic abundances, except for car- 
bon which is discussed below, and consider the reduced effec- 
tive cross-sections of atoms condensed in dust grains. The solar 
and cosmic abundance of carbon has long been considered to 
be around 400 ppm (e.g. Holweger 200 1) 5 , but Allende Prieto 
et al. (2002) re-determined the photospheric solar abundance 
to be 245 ppm, while Lodders (2003) suggested a protosolar 
abundance of 288 ± 30 ppm which is based on Allende Prieto's 
photospheric value and accounts for some elemental settling in 
the Sun. Measurements of the gaseous carbon abundance in the 
tenuous warm interstellar medium yield C/H 140 - 161 ppm 
(Cardelli et al. 1996; Sofia et al. 2004). There is a dispute 
about how much carbon is needed to explain the extinction 
properties of dust grains. Synthesis of the extinction curve for 
translucent lines of sight in terms of the properties of graphite, 

4 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec 

5 We quote abundances in parts per million (ppm) relative to hydro- 
gen. C/H=245 (400) ppm corresponds to a logarithmic abundance of 
8.39 (8.60). 
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amorphous carbonaceous grains, and silicate grains suggests 
that a minimum of 180ppm, and more typically 250 ppm, of 
carbon is needed to form grains with the observed properties 
(Clayton et al. 2003; Zubko et al. 2004), which suggests a car- 
bon abundance in the general interstellar medium somewhat 
larger than protosolar. The earlier proposition of an interstellar 
carbon abundance lower than solar photospheric (e.g. Wilms et 
al. 2000) has become unlikely with the reduction of the solar 
value. We generously condense the available information into a 
range of the total interstellar carbon abundance of 320 + 80 ppm 
of which ~ 60% is in dust grains. The size distribution of the 
grains is a subject of dispute, too, since the ones deduced from 
the extinction curves in different directions (e.g. Clayton et al. 
2003; Zubko et al. 2004) differ from that of grains of inter- 
stellar origin picked up in interplanetary space (Frisch et al. 
1999). Larger grains produce less extinction and will be less 
efficient X-ray absorbers. We give preference to the size distri- 
bution derived from extinction curves and use a MRN distribu- 
tion (Mathis et al. 1977; Clayton et al. 2003) with a cutoff at a 
maximum grain radius of 0.5 /itn as our best estimate to repro- 
duce the effects of carbonaceous grains. The influence of this 
parameter on our results is moderate. 

We fix the degrees of ionization of the interstellar hydro- 
gen and helium along the line of sight to our three stars. For 
HZ43 A, we use x H = 0.10 and x He = 0.40 (Wolff et al. 1999). 
For the other two stars, we adopt xh = 0.25 and xe e = 0.40, 
as determined from pickup ions of local interstellar origin 
and other studies (Slavin & Frisch 2002; Gloeckler & Geiss 
2004). Doubly ionized helium is disregarded following Slavin 
& Frisch (2002). 

4.3. The LETG+HRC spectrum of RXJ1856 

We followed Burwitz et al. (2003), Braje & Romani (2002), 
and Pons et al. (2002) in fitting the observed spectrum of 
RX J 1856 from the optical to the X-ray regime by the sum of 
two blackbodies with temperatures kr spot and kr star , of which 
the hotter one (kr spot ) describes the X-ray emitting spot. We 
provisionally interprete the cooler component (kr stal ) as the 
emission from the annulus representing the projected neutron 
star surface outside the spot 6 . 

We fit the summed model spectrum for the first to the 
sixth diffraction orders to the Chandra LETG+HRC-S count 
rate spectra of RX J 1856 in the positive and negative disper- 
sion directions. Fig. 4 shows the corresponding summed pho- 
ton spectrum, which allows to judge the strength of the inter- 
stellar carbon absorption edge more clearly than the count rate 
spectrum with its superposed effective area structure. The pho- 

6 Presently, the nature of the cool component is still disputed and so 
is the angular stellar radius derived from the blackbody assumption. 
A model in which an underlying blackbody flux is reprocessed in a 
thin surface layer of hydrogen can account for an enhancement of the 
optical flux relative to that at shorter wavelengths, but only by a factor 
of at most two (S. Dreizler, private communication) and, of course, 
it can not exceed the blackbody flux at its own temperature. A way 
to dissociate the optical component from the neutron star as a light 
source would be the presence of optical light emitting circumstellar 
matter. 
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Fig. 4. Chandra LETG+HRC spectrum of RXJ1856 binned to 0.5A, 
derived with the effective area at wavelengths longwards of the car- 
bon edge corrected to fit an abundance of carbon in the interstellar 
medium of a(C) = 245 - 400 ppm with 60% in dust grains. The two- 
temperature fit shown is for 320 ppm, with the upper level of the step 
indicated for 245 and 400 ppm, too. The fit has a^ 2 = 97.7 for 1 15 dof 
(see Sect. 4.5). The model parameters are quoted in Table 2. 



ton spectrum is created from the observed spectrum by sub- 
tracting the 2 nd to 6 th order model count rates from the data (in 
cts s _1 A _1 ) and dividing by the first-order effective area. The 
expected strength of the absorption edge is obtained if the ef- 
fective areas shortwards of 43 .7 A are accepted and longwards 
are moderately reduced. We perform fits to the wavelength re- 
gion 15-39 A and anchor the fit at long wavelengths where the 
higher orders dominate. This yields a data/model ratio differ- 
ent from unity in the 44. 0^-8. 5A range which we interprete 
as the mean factor a\ by which the effective area in this in- 
terval has to be adjusted in order to obtain a perfect fit with 
these wavelengths included. From best fits with a(C) between 
245 and 400 ppm and 60% of the carbon in dust grains, we ob- 
tain ai = (A) /(A ) = 0.875 ± 0.012, where the angular brack- 
ets refer to averages over the 44-48. 5A interval. An additional 
and larger uncertainty arises from the errors in the fit param- 
eters, notably kr star , determined for an increase in x 2 by +1 
over its minimum v 2 . =38.7 for 44 dof in the 15-39A range. 

^ mm c 

Merging the two results yields ce\ = 0.871 + 0.038. This fac- 
tor applies within its error to both the negative and the posi- 
tive dispersion directions. Figure 4 shows the derived summed 
spectrum of both sides using effective areas A = a(A) A where 
the wavelength-dependent adjustment factor a is derived be- 
low and is forced to match a\ between 44 to 48. 5 A. Ourphoton 
spectrum is very similar to that presented by Drake et al. (2002, 
see their Fig. 1), who evidently used similarly adjusted effec- 
tive areas. Our final fit with the adjusted areas extends over the 
spectral range 15-74 A and has x 2 = 97.7 for 1 15 dof. This 
spectral range avoids the shortest wavelengths with low flux 
and the longest wavelengths where the higher orders become 
important and finally dominate. 
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Fig. 5. Simultaneous fit of RX J1856, HZ43 A, and Sirius B in the wavelength ranges marked by vertical dotted lines (see Sect. 4.4.2). The 
LETG spectra binned to 0.5A are shown as data points, the corresponding best-fit models as solid curves, and the first-order contributions as 
dashed curves. The area correction function a is shown at the top. It converts the nominal LETG+HRC-S first-order effective area A of the 
November 2004 release into the adjusted area A used in this paper. Systematic uncertainties in a are indicated by error bars at 46, 70, 90, and 
125A. The steps in the count rate spectra of HZ43 A and RX J1856 at 49 and 69 A result from the detector gaps. Sirius B was observed off axis 
and its gaps are located differently (see text). 



4.4. Simultaneous fit to the spectra of several stars 
AAA. General approach 

The method for the determination of the area correction func- 
tion a(A) is based on the simultaneous fit to n > 2 stars with 
N data points each, centered at wavelengths A\ to A^. Because 
of a lack of suitable calibrators, we apply the method for n — 2 
only. The standard^ 2 which expresses the differences between 
the data with statistical errors cry and the model contribu- 
tions OTyk in orders k=l to 6 (all in cts s _1 A _1 ) is 



X 



EE 

j=i i=i 



(yij - Hk=l '«ijk) 2 



(1) 



with orders higher than k=6 entirely negligible for the objects 
considered. As usual, the model contributions are calculated as 



where f\(A\) is the incident spectral flux of star j in pho- 
tons cm~ 2 s _1 A~', £k(^i) the efficiency of the k-th order relative 



to the first order, and m ., are the model contributions calcu- 

ijk 

lated with the nominal effective areas. We correct the observed 
count rates for the respective model contributions in the second 
to sixth order to yield yi. - yij — 2^=2 a (Ai/k)m < L. Since a(/t;/k) 
is involved, the transformation requires, in principle, an itera- 
tive approach. While this is easily implemented in the iterative 
procedure to minimize x 2 , the influence of this procedure on 
the derivation of the c; = a(Ai) is minute because (i) the ad- 
justed areas affect the second and third orders only at A > 88A 
and A > 132 A, respectively, (ii) RX J 1856 is not affected be- 
cause we use its data only at A < 74A (see below), and (iii) the 
contributions of the higher orders to the spectra of HZ43 A and 
Sirius B are close to negligible. Hence the determination of the 
ct\ is straightforward. We substitute the yy in Eq. (1) by the yl , 
form the first-order count rate ratio of data vs. nominal model 

r ij = y[\l m \\ w ^ tn errors Pi] - an d find that the a\ are 

determined by minimizing the differences between the individ- 
ual ratios r« , 



X' 



=zz 

j=l 1=1 



0'', 



A) 2 



n N 

EE 

j=l 1=1 



Oij - ad 
p 2 . 



(3) 
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For n — 2 stars, we compute or; as the weighted mean of r ;i and 



an = 



p |ni +p|n2 



and obtain the final form of the^f 2 to be minimized as 

(ni-ra) 2 



(4) 



(5) 



i=l ^il 1 n2 

The expectation value of x 2 of Eq. (5) is N - p, since with 
2AT data points we determine p parameters in the models m\ 
and ni2 and N values of a\. The approach becomes insensitive 
to the a(A{) if the contributions of the higher orders at A\ be- 
come large, because the relative errors of the y'.. increase cor- 
respondingly. Minimizing x 2 with the a\ included may only be 
marginally stable if the parameters allow a bulging deforma- 
tion of both model spectra and a runaway of the a\. The white 
dwarf models are resistant against such deformation, however, 
and a(A\) is, furthermore, fixed to a\ = 0.87 in the interval 
44^-8. 5 A. In our application of the method, the fitted spectra 
of HZ43 A and RX J 1856 are tied down in the optical and at 
~ 46A, while the spectrum of Sirius B is tied to HZ43 A. 

We have opted to fit the sum of the model contributions 
for the positive and negative dispersion directions to the sum 
of the respective observed spectra. The so-defined model spec- 
trum allows us to derive the combined correction function a 
from the combined observed spectrum as well as the individual 
a-neg and a pos from the observed spectra for the negative and 
positive dispersion directions, respectively. The combined a is 
the area-weighted mean of a neg and ff pos , 



aA 



o 



,0 



(6) 



with A = Aj] eg + Ap OS . Multiplication of the nominal area with 
the respective a yields the adjusted areas A neg 
Apos — ^posAp and A = oA^ . 



a negA„ eg , 



4.4.2. Combined fit to RXJ1856, HZ43A, and SiriusB 

The number of stars with sufficiently well-known incident 
spectra is extremely restricted. We fit, therefore, RXJ1856 
alone for 15 < A < 39 A, RXJ1856 and HZ43A simultane- 
ously for 44 < A < 74A, and Sirius B and HZ43 A simultane- 
ously for 74 < A < 160A. We calculated a grand x 2 f° r tne 
combined fit of the three models to the three observed spectra 
with 518 data points including three points in the visual and 
the two 485 A fluxes of HZ43 A and Sirius B (the latter chosen 
to reproduce the Hebrard et al. (1999) value of Nm within its 
error). The number of free parameters includes the 12 listed 
in Table 2 plus the 232 values a\. Hence, there are 274 dof. 
Figure 5 displays the results of the simultaneous fit to the three 
stars and Table 2 gives the best-fit parameters. The quoted 1-cr 
errors are derived from variation of the respective parameter 
with all other parameters kept free and refer to Ax 2 = +1. The 
only parameter kept fixed in the fit was the gravity of Sirius B, 
logg = 8.6, which is reasonably well determined from the an- 
gular radius, the Hipparcos distance and the gravitational red- 
shift (Barstow et al. 2005; Holberg et al. 1998). The grand x 2 



Table 2. Parameters of HZ43 A, SiriusB, and RXJ1856 based on the 
simultaneous fit of our model spectra to the LETG+HRC count rate 
spectra in the wavelength intervals given. The quoted 1-cr (Ax 2 = +1) 
errors are correlated and derived from fits with the other parameters 
for each object kept free. The letter / indicates: fixed. 



Parameter 



Value±Error 



(a) HZ43 A (A = 45 - 160A) 
r eff (K) 

R 2 /cP (10- 23 ) 



N m (10 17 cnr 2 ) 



(b) SiriusB (A 
Ta (K) 
logg 

R 2 /d 2 (10- 21 ) 
A H i(10 17 cnr 2 ) 

(c) RXJ1856 (A = 15 
kr spot (eV) 
kr star (eV) 
RJd (km/pc) 
R 2 /d (km/pc) 
N m (10 20 cnT 2 ) 



74- 160A) 



-74A) 



51126 ±660 

7.90 ± 0.08 
3.011 ±0.038 

8.91 ± 0.37 

24923 ±115 

8.6 / 1 
4.877 ± 0.036 

6.5 ±2.0 2 

62.83 ± 0.41 
32.26 ± 0.72 
0.0378 ± 0.0003 
0.1371 ±0.0010 
1.10 ±0.03 



1 Based on Barstow et al. (2005); Holberg et al. (1998) 

2 Hebrard et al. (1999). Our fit is required to stay within the 1-cr error. 



for the simultaneous fit is 282.9 for 274 dof (xl = 1.03) which 
is close to perfect. 

The self-consistently determined multiplicative correction 
function a(A() = A/A () for the sum of negative and positive dis- 
persion orders is shown at the top of Fig. 5. It converts the nom- 
inal LETG+HRC-S effective area A of the November 2004 re- 
lease to the area A used in this paper and is valid for the best-fit 
log g=7.90 of HZ43 A. It has been smoothed with a box car 
over three wavelength bins (i.e., a total width of 1.5 A) in order 
to reduce the influence of statistical fluctuations in the individ- 
ual spectra. The correction is consistent with unity at A < 44A, 
but reaches a low of about 0.77 at 65 A, consistent with Pease's 
et al. (2003) expectation that the error in the 60-80A range 
might reach 25%. At longer wavelengths, a(A\) moves closer 
to unity again. For logg values of HZ43 A other than 7.90, 
the whole curve moves up and down a bit with a pivot point 
at a(46A) = 0.87 ± 0.04. In order to mark the extent of this 
variation, we have included the systematic errors at 46, 70, 90, 
and 125 A with an upper (lower) cross bar signifying the lo- 
cation of a(Ai) for a, = 0.91(0.83) and logg=8.00 (7.80) for 
HZ43 A, respectively. None of the other parameters has a sim- 
ilarly large influence. In summary, we conclude that the nec- 
essary adjustment to the LETG+HRC-S effective area relative 
to that at A ^ 40A is pronounced in the 50-80A region, but 
small longwards of 85 A. Its systematic uncertainty is +5% at 
A <: 80A from the uncertainty in logg of HZ43 A and another 
5% at all wavelengths A <: 44A from the error in a\. The fact 
that log g of Sirius B is better defined than that of HZ43 A does 
not help because the spectrum of Sirius B does not extend to 
44A where it could be tied down. 
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Fig. 6. Photon spectra of RXJ1856, HZ43 A, and Sirius B based on 
our analysis of the LETG+HRC-S observations (data points) and best 
fit model spectra with parameters as given in Table 2 (solid curves). 
The spectrum of RX J 1 856 is the same as in Fig. 4 with the two black- 
body components shown separately as dashed curves. Longwards of 
78A, we have included the EUVE spectra of the two white dwarfs 
(green/gray curves). 

Our derived areas for the negative and positive dispersion 
directions, A neg = a n eg^eg and ^pos = a P o S Ap OS , are close to the 
SRON areas suggested previously by J. Kaastra (Pease et al. 
2000, and J. Kaastra, private communication) between 44 and 
80A and close to the nominal CXC areas at longer wavelengths. 
Our results are provided in Tables Al to A3 7 for the preferred 
log g of HZ43 A of 7.90 as well as for 7.80 and 8.00, indicating 
the systematic errors of our fit. The following systematic uncer- 
tainties are not included in these Tables: (i) the normalization 
uncertainty at ~ 46A as expressed by a\ = 0.87 ± 0.04; (ii) the 
possible degradation in the HRC response with time; and (iii) 
the absolute error in the nominal areas A at short wavelengths. 

As noted above, our empirically adjusted effective areas 
are strictly valid only for the epoch JD 2452 130. The possi- 
ble variation of the HRC sensitivity since the begin of the mis- 
sion is small, however, 0.015 in a{A\) at long wavelengths and 
possibly as much as 0.050 shortwards of about 50A (Fig. 1). 
Nevertheless, the possibility of a slow degradation in A neg and 
A pos should be kept in mind: it ultimately limits the accuracy 
of parameters determined from fits to the long wavelength sec- 
tions of LETG+HRC spectra. 

4.4.3. The photon spectra of HZ43A, Sirius B, and 
RXJ1856 

The incident photon spectra of the three stars fitted simultane- 
ously are obtained by dividing the count rate spectra of Fig. 5 
by our effective areas. Their errors are dominated by the sys- 

7 Tables Al to A3 are available only electronically at 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/. 



Table 3. Soft X-ray fluxes as defined by the best-fit absorbed photon 
spectra in Fig. 6. The numbers in brackets are the systematic errors of 
the LETG+HRC-S area correction in percent (see text). 
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2.14 x 10~ 3 
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2.49 x 10- 4 


(6) 
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2.860 x 10 
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2.96 x 10~ 19 3) 
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2.558 x 10 


-14 









" The Compton effect may depress this flux by ~ 1%. 

2) The unabsorbed flux is 6.27 x 10~ 17 ergcm _2 s _1 A _1 

3) The unabsorbed flux is 3.18 x 10 - ' 9 ergcm _2 s _1 A _1 



tematic errors as discussed in the last Section. Fig. 6 shows 
the resulting photon spectra with the systematic uncertainties 
indicated again by additional error bars at 60 and 125 A. The 
spectral fluxes at selected wavelengths are listed in Table 3 to- 
gether with their systematic errors arising from a in percent 
given in brackets. The complete spectra are provided in Tables 
A4 to A6 8 . The systematic errors listed near the end of the last 
Section under (ii) and (iii) are not included. The bottom three 
lines of Table 3 contains the optical fluxes to which the spec- 
tra were normalized. The fourth line from the bottom gives the 
predicted fluxes at 1300A. In the case of HZ43 A, the quoted 
flux agrees to better than 0.3% with that of the HST/STIS spec- 
trum. For Sirius B, the 1300A model flux agrees with that pre- 
dicted by Vennes & Dupuis (2002) and both exceed the IUE 
flux which is known to be too low (Massa & Fitzpatrick 2000). 

4.5. Results for individual stars 

We now discuss the results for the individual stars taking the 
systematic errors as far as possible into account. These fits are 
characterized by x 2 =97. 7 for RX J1856 in the interval 15-74A 
with 115 dof, * 2 =183.9 for HZ43 A in 45.5-160A with 229 
dof, and^ 2 =209.4 for SiriusB in 61.5-160A with 197 dof. All 
fits are substantially improved over those without the area cor- 
rection and can all be considered 'good' in a statistical sense. 
This result is not subject to the binning chosen as in previous 

8 Tables A4 to A6 are available only electronically at 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/. 



10 



K. Beuermann et al.: Soft X-ray standard stars HZ43 A, Sirius B, and RX J185635-3754 



1.1 



1.0 



0.9 



C-edge 



HZ43A 



mm 



o.e 



40 



60 

Wavelength (A) 



100 



Fig. 7. Data vs. model flux ratio of the spectrum of HZ43 A for A > 
43. 5A. The maximum depression of the spectrum by the Compton 
effect near the carbon K-edge at 43. 7A is at the percent level. 



fits to the LETG+HRC spectrum of RXJ1856 (e.g. Braje & 
Romani 2002). The combined^ 2 for HZ43A and Sirius B is 
393.3 for 426 dof and the slight disparity between the two stars 
results form the better statistics of HZ43 A which drives the fit. 



4.5.1. HZ43A 

The errors of the individual parameters listed in Table 2 are, 
in part, highly correlated. This is particularly true of T e ff and 
log g for HZ43 A. The corresponding error ellipses for the 
68%, 90% and 95% confidence limits are overplotted in Fig. 3. 
Gratifyingly, the best-fit values are in good agreement with re- 
sults obtained from the optical and ultraviolet regimes. The 
large extent along the major axis of the ellipses reflects the 
correlation between the results for T e s and log g. The soft X- 
ray flux stays essentially constant along the major axis and 
the narrow width in the perpendicular direction signifies the 
rapid deterioration of the fit if the model does not meet the 
observed flux. Hence, the soft X-ray flux in Table 3 is better 
defined by the fit than either T e g or log g. Parameter combina- 
tions like those in the lower left of Fig. 3 can safely be excluded 
since they yield much too low a flux, parameter combinations 
in the upper right would yield too large a flux. This is an impor- 
tant result for the future calibration of instruments sensitive in 
this spectral regime: it is more appropriate to quote soft X-ray 
fluxes than an effective temperature. 

The angular radius R/d of HZ43 A fixes the distance d for 
a given radius R. Adopting Wood (1995) stellar models with a 
thick hydrogen envelope provides an independent relation and 
yields M = 0.64 + 0.05 M @ , R = (1.04 + 0.08) x 10 9 cm, and 
d = 61.3 + 5.5 pc with errors which account for the systematic 
uncertainty due to the full gravity range logg = 7.90 + 0.10. 
The gravitational redshift predicted for these parameters is 
v g = 28 ± 4 km s -1 . The observed i/ g = 30.1 ± 1 .5 km s" 1 (Reid 
1996) suggests the slightly more restricted range M = 0.67 + 
0.02 M Q (Dupuis et al. 1998), with R = (9.92 + 0.34) x 10 8 cm, 
logg = 7.95 + 0.04, and d = 59 + 2pc (1-cr errors). Both 
distance values are consistent with the trigonometric parallax 
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160 



Fig. 8. Ratio of the soft X-ray fluxes in HZ43 A and Sirius B based 
on the Chandra LETG+HRC spectra analyzed here and the EUVE 
archive spectra listed in Table 1. The results are valid for the cali- 
bration using log g(HZ43 A)=7.90 and would shift up or down by the 
systematic error shown at 128 A for logg lower or higher by 0.1, re- 
spectively. 



Tabs = 15.3 + 2.9 mas (d = 65+jgpc) from the Yale Catalog of 
Trigonometric Parallaxes (van Altena et al. 2001). 

In principle, the LETG+HRC data of HZ43 A can help 
to decide whether the short-wavelength spectra of hot white 
dwarfs are steepened by the fact that electron scattering is no 
longer in the Thomson limit and must be correctly described 
by the Compton effect which degrades the energies of the scat- 
tered photons. In Fig. 7, we show the data/model ratio of the 
photon spectrum of HZ43 A from Fig. 6 for A > 44A. The 
data points are all more or less consistent with a ratio equal to 
unity: 1.0005 + 0.0014 for A = 65 - 100A; 0.995 ± 0.011 for 
A = 49- 58A; and 0.959 + 0.027 for A = 44 - 49A. The glitch 
in the ratio near 60A may be a relic of differences between the 
positive and negative dispersion directions. We can not exclude 
an effect at the percent level given the remaining systematic er- 
rors, but the effect may as well be still smaller as suggested by 
Suleimanov et al. (2006). 

4.5.2. Sirius B 

The angular radius for Sirius B in Table 2 falls between the val- 
ues quoted by Barstow et al. (2005) for the HST/STIS G430L 
and G750L spectra. At the Hipparcos distance of 2.637 + 
0.01 lpc, the radius of the white dwarfism = (5.683 + 0.019) x 
10 s cm = (8.165 + 0.028) x 10 3 R Q , where the error includes 
the systematic error in the LETG+HRC-S area correction and 
is insensitive to logg of Sirius B itself. Combined with the 
gravitational redshift of Barstow et al. (2005), v gr = 80.42 + 
4.83 km s~', we obtain amass of M = 1.032 + 0.062 M @ , con- 
sistent with the astrometric solution of Gatewood & Gatewood 
(1978). The implied gravity logg = 8.627 + 0.026 slightly ex- 
ceeds the spectroscopic value of Barstow et al. (2005), but stays 
within Holberg's et al. (1998) 1-cr error. Conversely, the effec- 
tive temperature of 24923 + 1 1 5 K (systematic error) falls below 
Barstow's spectroscopic value of 25193 + 37 K (statistical er- 
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ror) derived from a Balmer line fit to the HST/STIS spectrum. 
The possible systematic error quoted by Barstow et al. (2005) 
is 350K and encompasses our result 9 . A spectrum calculated 
with our code for Barstow's T e ff and \ogg adjusted to the same 
visual HST/STIS flux as above yields a predicted LETG+HRC 
spectral flux 14% higher than our spectrum in Fig. 6, about 
twice its log ^-derived systematic error of 7%. 

4.5.3. RXJ1856 

The high quality (xl = 0.87) of the two-blackbody fit to the 
LETG+HRC spectrum of RX J 1856 in Figs. 4 and 6 demon- 
strates that the deviation from this model is minute. The model 
is normalized to the extinction-corrected optical flux of van 
Kerkwijk & Kulkarni (2001) and provides an almost perfect 
fit to all available data. Blackbody temperature and angular 
radius of the hotter component, kr spot = 62.8 + 0.4 eV and 
0.0378 km pc -1 (Table 2), agree well with previous determi- 
nations (Burwitz et al. 2001, 2003; Braje & Romani 2002; 
Drake et al. 2002; Pons et al. 2002). The temperature kr star = 
32.3 + 0.7 eV of the cooler component is determined here more 
accurately than before as a result of our elimination of some of 
the deficiencies in the LETGS effective area calibration at long 
wavelengths. The angular radius of this component (Table 2) 
refers to the outer radius of a blackbody-emitting annulus sur- 
rounding the hot spot. The quoted errors account for the sys- 
tematic uncertainty in the area correction and indicate that the 
main error source for the radii still lies in the distance uncer- 
tainty. The excellence of the two-blackbody fit does not imply 
that RXJ1856 actually radiates as a blackbody (see, e.g., the 
discussion in Kargaltsev et al. 2005). E.g., a graybody with its 
different limb darkening law would yield a different radius than 
the disk of uniform brightness. However, whatever the emission 
properties, the spectral flux integrated over the surface must 
nevertheless agree with that of the blackbody model within 
the statistical uncertainties of the X-ray and optical/ultraviolet 
data. Hence, the blackbody fit is, in fact, quite restrictive. For 
the purpose of using RXJ1856 as a soft X-ray calibrator the 
blackbody assumption is entirely satisfactory and need be re- 
placed only if future observations reveal fine structure in the 
spectra not accounted for by this approximation. At present, 
we agree with (Braje & Romani 2002) that there is no evi- 
dence for a more complex temperature distribution beyond the 
dichotomy of base and spot temperature. 

4.6. Calibration of the EUVE SW spectrometer 

It is instructive to compare our photon spectra of HZ43 A and 
Sirius B with the short-wavelength spectra of these stars in 
the EUVE archive (Table 1). The mean EUVE spectra are in- 

9 Note that our systematic error of 1 15 K does not include the uncer- 
tainty in the visual flux of Sirius B: a +0.01 mag change (1% decrease 
in the adopted 4600A flux) raises our value of the effective tempera- 
ture of Sirius B by 14 K or 0.057%. It is in the nature of our simul- 
taneous fit that at the same time the effective temperature of HZ43 A 
is slightly reduced (by 0.050%). The soft X-ray fluxes stay practically 
unchanged. 
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Fig. 9. ROSAT PSPC spectra of HZ43 A and Sirius B in either the open 
configuration or with boron filter (B). The abscissa is in channels with 
one channel corresponding to an apparent photon energy of 10 eV. The 
model spectra based on the fits to the LETGS observations are shown 
as solid curves, the PSPC observations as data points with statistical 
errors. HZ43A possesses a bremsstrahlung component which com- 
mences at channel 65 and is included in the fit. 



eluded in Fig. 6 spectra and are seen to fall below the Chandra 
LETG+HRC spectra in our calibration. Fig. 8 shows the ratios 
of Chandra vs. EUVE fluxes. They average 1.14 and 1.16 for 
HZ43 A and Sirius B, respectively. The systematic error of 7% 
from the uncertainty in log g of HZ43 A would shift both num- 
bers up or down. This comparison suggests that the effective 
area used for the short-wavelength section of the EUVE archive 
spectra is too high by a factor of 1 . 15 ± 0.07 (systematic error). 
A similar discrepancy between their best-fit models and the 
EUVE short-wavelength fluxes was noted previously by Sing 
et al. (2002). Interestingly, they found agreement of the mod- 
els with the EUVE fluxes at longer wavelengths, a result which 
has implicitely been adopted in our approach of determining 
the hydrogen column density by fitting the model spectrum to 
the EUVE long-wavelength flux at 485 A. Our factor of 1.15 at 
the short wavelengths happens to be close to that by which the 
EUVE short-wavelength effective area determined on ground 
was corrected upward after an in-flight calibration using theo- 
retical spectra of the white dwarfs HZ43 A, GD71, and GDI 53 
(Boyd et al. 1994). In retrospect, that correction seems to have 
been unnecessary. 

4.7. ROSAT PSPC calibration 

We test the ROSAT PSPC calibration by folding the photon 
spectra of Fig. 6 through the PSPC detector response and 
compare them with the observed low-resolution PSPC spectra. 
We consider here the spectra of HZ43 A and Sirius B which 
were observed early in the mission (PV and AO-1 phases) and 
about two years into the mission (AO-2), before and after the 
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Table 4. Observed PSPC count rates in channels 1 1-60 and LETGS- 
predicted count rates based on the photon spectra in Fig. 6. 'O' and 
'B' in column 3 refer to observations in the open configuration and 
with the boron filter, respectively. 











Count rate 


Object 


Date 


Filter 


Exp. 


observed 


predicted 




Start 




(ks) 


(cts/s) 


(cts/s) 


HZ43A 


900621 


O 


21.5 


66.08 ± 0.30 


65.38 ± 2.70 




910619 


B 


21.6 


8.04 + 0.03 


7.35+0.28 


Sirius B 


910315 


B 


3.3 


1.43 + 0.03 


1.38 + 0.05 


RXJ1856 


921111 


O 


6.3 


3.51+0.03 


3.33 + 0.17 



change in high voltage which affected the detector response 
(see Sect. 2.3 and Table 1). We have selected two long obser- 
vations of HZ43 A, one in the open configuration during the PV 
phase with PSPC-C and one with the boron filter in AOl with 
PSPC-B. Sirius B was observed only once with the boron filter 
in AO-1. We used the appropriate standard detector matrices 
DPvMPSPC-AOl and DRMPSPC and applied no corrections 
for gain shifts (see Sect. 2.3). The resulting PSPC spectra are 
shown in Fig. 9 and reproduce the observations exceedingly 
well. The relevant integrated count rates for energy channels 
1 1-60 are listed in Table 4. The errors given for the predicted 
count rates refer to the systematic uncertainty in our photon 
spectra and LETG+HRC-S area correction discussed above. 
The observed count rate of HZ43 A in the open configuration 
agrees perfectly with the prediction. Hence, we can safely con- 
clude that the calibration of the PSPC is correct and the sugges- 
tions of a miscalibration by Napiwotzki et al. (1993), Jordan et 
al. (1994), and Wolff et al. (1995, 1996) are unfounded. Close 
inspection reveals some small discrepancies between predic- 
tion and observation which find ready explanations in a variety 
of small necessary corrections. In the uppermost spectrum in 
Fig. 9, the HZ43 observation in the open mode, the observed 
spectrum is slightly shifted to lower energy channels. This ef- 
fect is caused by the gain depression by some 4% in the small 
central focal spot (Snowden et al. 2001). The observational 
'wobble' mode moves the source forth and back across this 
central depression and causes a mean gain reduction of about 
1.5% which accounts for the observed shift. Of a different ori- 
gin is the slight excess in count rate for the observations with 
boron filter (lower two spectra in Fig. 9). A probable explana- 
tion is provided by an actual thickness of the boron window 
about 4% lower than nominal. 

The remaining small calibration uncertainties of the PSPC 
and the calibration of the high resolution imager (HRI) will be 
addressed elsewhere. All these effects are at the few percent 
level and provide no basis for the allegation of a miscalibra- 
tion. In summary, the internal consistency of our results sug- 
gests that the PSPC yields a correct measurement of the soft 
X-ray flux and the ground calibration of the detector at soft X- 
ray energies is reliable to better than 10%, as concluded also by 
Snowden et al. (1995). Given the large number of observations 
of all types of celestial sources, both in the All-Sky Survey and 
in pointed observations, we judge this a reassuring statement 



of general interest. One can only speculate on the basis of the 
reports about a major miscalibration of the PSPC. The use of 
model spectra which provided perfect fits to observations in the 
optical and ultraviolet but lacked flux in the soft X-ray regime 
is a possibility. 

In passing, we confirm that HZ43A possesses a 
bremsstrahlung component with a temperature of 0.6 keV and 
an emission measure of 3 x 10 52 cirT 3 , which is seen to com- 
mence near channel 65 in Fig. 9. It was previously reported by 
O'Dwyer et al. (2003) and probably originates from the sec- 
ondary star. We have checked the zeroth order images of all 
individual LETGS observations for a signal at the position of 
the secondary (separation 3 arcsec, position angle 280°), but the 
point spread function of the white dwarf image still contributes 
a flux about an order of magnitude larger than expected for the 
secondary even if the position of the secondary is favorably lo- 
cated between the diffraction spikes. Hence, an origin of the 
hard X-rays from the secondary is plausible but remains un- 
proven. 

5. Conclusions 

We have established the neutron star RXJ1856 and the hot 
white dwarfs HZ43 A and Sirius B as standard stars in the soft 
X-ray regime by fitting the LETG+HRC spectra of all three 
stars simultaneously with the best available model spectra for 
the white dwarfs and a two-blackbody model for the neutron 
star. Our method ties the soft X-ray spectra of the two white 
dwarfs to that of the neutron star RX J 1856 and thereby to be 
better calibrated effective areas of the Chandra LETG+HRC 
spectrometer shortwards of 40A. The spectrum of HZ43 A (and 
indirectly that of Sirius B) is thus fixed at soft X-ray and at op- 
tical wavelengths, a procedure which avoids the wagging-tail 
problem of previous calibration attempts using white dwarfs. 
Our second major result is the self-consistently determined cor- 
rection to the LETG+HRC-S effective areas longwards of 44A. 
We find that the areas for both dispersion directions need a 
maximum correction of about -23% at 65 A and are correct 
to within ~ 10% at wavelengths exceeding 80A. These correc- 
tions are valid relative to the short-wavelength effective areas. 
The largest remaining uncertainty, which is not addressed here, 
lies in the absolute calibration of the LETG+HRC-S at short 
wavelengths. Our photospheric parameters of the three stars 
are in good agreement with those derived from independent 
optical/ultraviolet and soft X-ray studies and, in part, improve 
on them. The soft X-ray flux of HZ43 A can be produced by 
a range of T e g-\ogg combinations and is better defined than 
either of these two parameters. 

We find that the calibration of the EUVE short-wavelength 
spectrometer differs from that of the Chandra LETG+HRC-S 
by about 15 + 7% in the sense that the original ground cal- 
ibration seems to have been more adequate than the adopted 
in-flight calibration. 

We demonstrate that the ground calibration of the ROSAT 
PSPC is correct to within a few percent and that previous re- 
ports of a major miscalibration are unfounded. The observed 
PSPC spectra of our three stars agree within a few percent with 
the LETGS-based soft X-ray spectra folded through the PSPC 



K. Beuermann et al.: Soft X-ray standard stars HZ43 A, Sirius B, and RX J185635-3754 



13 



response. This internal consistency, in turn, supports our ad- 
justment of the LETG+HRC-S calibration at long wavelengths 
and the basic correctness of the nominal calibration at short 
wavelengths. 
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